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ABSTRACT 

Swift/BAT has detected ~ 200 long-duration GRBs, with redshift measurements for ~ 50 of them. We 
derive the luminosity function (<I> HL ) and the local event rate (p^ L ) of the conventional high luminosity 
(HL) GRBs by using the z-known Swift GRBs. Our results are generally consistent with that derived 
from the CGi? 0/BATSE data. However, the fact that Swift detected a low luminosity (LL) GRB, GRB 
060218, at z = 0.033 within ~ 2 year of operation, together with the previous detection of the nearby 
GRB 980425, suggests a much higher local rate for these LL-GRBs. We explore the possibility that 
LL-GRBs as a distinct GRB population from the HL-GRBs. We find that pfr L is ~ 325±f% Gpc" 3 yr"\ 
which is much higher than p$ h (1.12± j% Gpc~ 3 yr" 1 ). This rate is - 0.7% of the local Type Ib/c SNe. 
Our results, together with the finding that less than 10% of Type Ib/c SNe are associated with off-beam 
GRBs, suggest that LL-GRBs have a beaming factor typically less than 14, or a jet angle typically wider 
than 31°. The high local GRB rate, the small beaming factor, and low luminosity make the LL-GRBs 
distinct from the HL-GRBs. Although the current data could not fully rule out the possibility that both 
HL- and LL-GRBs are the same population, our results suggest that LL-GRBs are likely a unique GRB 
population and the observed low redshift GRB sample is dominated by the LL-GRBs. 

Subject headings: gamma-rays: bursts — gamma-ray: observations — methods: statistical 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) and supernovae (SNe) are 
two of the most violent explosions in the Universe. The 
connection between long duration GRBs and SNe was pre- 
dicted theoretically (Colgate 1974; Woosley 1993), and has 
been verified observationally through detecting spectro- 
scopic features of the underlying SNe in GRB 980425/SN 
1998bw (Galama et al. 1998; Kulkarni et al. 1998), GRB 
030329/SN 2003dh (Stanek et al. 2003; Hjorth et al. 

2003) , GRB031203/SN 20031w (Malesani et al. 2004), and 
GRB 060218/SN 2006aj (Modjaz et al. 2006; Pian et al. 
2006; Sollerman et al. 2006; Mirabal et al. 2006; Cobb 
et al. 2006). In some other cases, red SNe bumps were 
observed in the late optical afterglow light curves (Bloom 
et al. 1999, 2002; Delia Valle et al. 2003; Fynbo 2004; 
see a comprehensive sample and references in Zeh et al. 

2004) . It is long speculated that long GRBs are associ- 
ated with deaths of massive stars, and hence, SNe (for 
recent reviews, see Zhang & Meszaros 2004; Piran 2005; 
Meszaros 2006; Woosley & Bloom 2006). This specula- 
tion was broken down by the recent observation of GRB 
060614, which is a long, nearby GRB without an accom- 
panied SN (Gehrels et al. 2006, Fynbo et. al. 2006, Delia 
Valle et al. 2006, Gal- Yam et al. 2006). Based on its close 
analogy with the short duration GRB 050724, Zhang et 
al. (2007a) argued that this event might be also powered 
by a merger of compact stars. They further suggest that 
the conventional long vs. short GRB classification scheme 
may be modified as Type I (from mergers) vs. Type II 
(from collapsars) GRBs, and that GRB 060614 belongs to 
Type I (see also Zhang 2006). 

Some authors have attempted to determine the lumi- 



nosity function ($) and the local rate (po) of long GRBs 
through fitting the log A — logP or V/U max distributions 
observed by CGRO/BATSE (Schmidt 2001; Stern et al. 
2002; Lloyd-Ronning et al. 2002; Norris 2002; Guetta et 
al. 2005) or through simulations (Lloyd-Ronning et al. 
2004; Dai & Zhang 2005; Daigne et al. 2006). They con- 
sider only high luminosity GRBs (HL-GRBs, with lumi- 

, and generally characterize <E> 



1 Gpc 3 yr 1 



nosity L > 10 9 erg s 
with a broken power law and obtain p^ L 
(e.g. Schmidt 2001; Guetta et al. 2004, 2005). Guetta 
et al. (2004) suggested that by extrapolating the LF of 
HL GRBs to low-luminosities with a broken power law, 
one gets Pq L as ~ 10 Gpc -3 yr -1 for low redshift GRBs 
atz< 0.17 (GRBs 980425, 031203, and 030329). How- 
ever, the fact that Swift detected GRB(XRF) 060218 at 
z = 0.033 (Mirabal et al 2006) within ~ 2 years of opera- 
tion, together with the previous discovery of GRB 980425 
at z = 0.0085 (Tinney et al. 1998) by BeppoSAX, suggests 
that the local rate of the GRB 060218-like low luminos- 
ity (LL) GRBs (p% L ) is 100 ~ 1000 Gpc~ 3 yr" 1 (Cobb 
et al. 2006; Pian et al. 2006; Soderberg et al. 2006b), 
much higher than p$ h (Schmidt 2001; Stern et al. 2002; 
Lloyd-Ronning et al. 2002; Norris 2002; Guetta et al. 
2004, 2005). This poses a great puzzle about the nature 
of of these LL-GRBs. Are these LL-GRBs from the same 
population as HL-GRBs (Nakamura 1999; Ioka & Naka- 
mura 2001; Yamazaki et al. 2003; Guetta et al. 2004) 
or from a sub-energetic GRB population (Kulkarni et al. 
1998; Soderber 2004b; Mazzali et al. 2006; Toma et al. 
2006) with a much higher rate than that of the conven- 
tional HL-GRBs? In this paper we attempt to address 
these questions through careful statistical analyses with 
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the redshift-known Swift GRB sample. Throughout the 
paper H = 71 km s _1 Mpc -1 , m = 0.3, and A = 0.7 
are adopted. 

2. DATA AND GRB SAMPLE 

The current GRB sample with rcdshift measurements 
contains more than 80 GRBs. Since both the observed 
luminosity and rcdshift distributions are instrument- 
dependent (Jakobsson et al. 2006; Zhang et al. 2007b), 
we use only the long GRBs detected by Swift to form 
a homogeneous sample. Swift has detected ~ 200 long 
bursts during the first 2 years of operation (from Decem- 
ber 2004 to November 2006). Among them 47 bursts have 
rcdshift measurements. We have excluded GRBs 060505 
and 060614 from our sample since they are suggested to 
be classified as a new type of GRBs (Gehrels et al. 2006; 
Gal- Yam et al. 2006; Fynbo et al. 2006; Delia Valle et 
al. 2006) or to belong to the short GRB group (Zhang et 
al. 2007a; Zhang 2006). We therefore obtain a Swift HL- 
GRB sample of 45 bursts. We collect the peak flux and 
spectral parameters of these GRBs from the literature or 
GCN Circular reports. It is well known that most broad 
band GRB spectra are well fitted by the Band function 
(Band et al. 1993). However, ~ 80% of the BAT spec- 
tra are fitted by a simple power law (F v oc v~ T ) (e.g., 
Zhang et al. 2007b). This effect is due the narrowness 
of the BAT band (15-150 kcV) and the faintness of the 
bursts (so that there are not enough counts to perform a 
Band-spectrum fit. Therefore, the BAT spectra are not 
the true spectra of the Swift GRBs. In order to correct 
the observed luminosity to a broad band (1 — 10 4 keV in 
this analysis) we have to estimate the parameters of the 
Band function for each GRB. Zhang et al. (2007b) devel- 
oped a new method to derive the parameters of the Band 
function by fitting the observational data with incorpora- 
tion of the observed spectral hardness ratio. They found 
a strong correlation between T PL derived from the BAT 
data and E p of the vf v spectrum (see also Zhang et al. 
2007a). Sakamoto et al. (2006) independently derived 
this correlation. We use this method to estimate the pa- 
rameters of the Band function for the GRBs in our sample, 
and then correct the observed luminosity to 1 — 10 4 keV 
band in the burst rest frame with the fc-correction pre- 
sented by Bloom ct al. (2001). The HL-GRB sample has 
a moderate size, with both a broad luminosity distribution 
[L lso = 2 x 10 50 ~ 4.6 x 10 53 erg s" 1 ) and a broad redshift 
distribution (z = 0.55 - 6.29). 

3. ANALYSIS METHOD 

The number of GRBs per unit time at redshift z <~ z+dz 
with luminosity L ~ L + dL is given by 
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where Rgrb(z) is the event rate (in units of Gpc~ 3 yr _1 ) 
as a function of z, (1 + z)~ 4 accounts for the cosmological 
time dilation, and 



dV 



is the comoving volume element at redshift z in a flat A 
CDM universe. We assume that -Rqrb follows the star 
formation rate as a function of redshift, and the parame- 
terized star formation model SF2 presented by Porciani & 
Madau (2001) is used 1 , 
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where po is the underlying local GRB rate per unit volume 
at z ~ 0, i.e., po = -RgrbU^o- We characterize $(£) as 



(4) 



where Lb is the break luminosity and $o is a normaliza- 
tion constant to assure the integral over the luminosity 
function being equal to unity. Considering an instrument 
having a flux threshold F t h and an average solid angle O 
for the aperture flux, the number of the detected GRBs 
after an observational period of T should be 



N = — I <f>(L)dL 



Rgrb(z) dV(z) 



l + z 



dz 



dz, (5) 
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where L max and L m ; n are taken as 10 54 and 10 45 erg s _1 , 
respectively, and z max for a given burst with luminosity 
L is determined by the instrumental flux threshold F t h 
through F th = L / An D 2 L (z max ) . 

We derive $(L; a\, ct2, Lt>) with the Swift GRB sam- 
ple by using the following procedure. First, we constrain 
ai, «2, and Lb by comparing the luminosity and redshift 
distributions predicted by the model to the observations 
(1-D Criterion). We calculate the L and z distributions 
for a given set of the parameters, and then measure the 
consistency between model predictions and observations 
with the K-S test (Press et al. 1999). The significance 
level of the consistency is described by a probability of 
Pk-s- A larger value of pk-s suggests a better con- 
sistency (the higher significance level). The consistency 
for both L and z distributions is therefore measured by 
p 1 K_s = Pk-s x Pk-S- Second, we derive po for a given $ 
by using N — N obs (A^-Criterion). In this analysis we con- 
sider only the GRBs detected by Swift/BAT. The aperture 
solid angle of BAT is 1.33, and the detection numbers of 
both HL- and LL-GRBs by Swift/BAT in 2 operation years 
is Nfil = 200 and N°l s = 1. The sensitivity curve of BAT 
in the 50-150 keV band (Band 2003) is adopted. Third, 
we use a 2-dimensional GRB number distributions in the 
[logL, log(z)]-plane to further examine the parameters 
and evaluate the influence of the observational biases of 
the sample (2-D Criterion). The 1-D Criterion constrains 
the parameters with the clustered bursts in the sample. 
The 2-D Criterion, instead, present a self-consistency test 
with the scattered data points by assuming that the burst 
distribution in the [logi, log(z)]-plane predicted by a rea- 
sonable <& should cover the observational data points at a 
3cr confidence level. 

4. ARE HL-GRBS AND LL-GRBS THE SAME POPULATION? 



dz #o (1 + 2)2 [0^(1 + 2)3 + 0^1/2 

1 We have considered the other two star formation models, SF1 and SF3. Our results are not sensitive to the selection of these star formation 
models since these models are almost the same at low rcdshifts. 
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We first study the constraints on $ HL with the the cur- 
rent z— known HL-GRB sample, and then analyze whether 
or not it can predict the detection rate of LL-GRBs by ex- 
trapolating $ HL down to a luminosity as that of GRBs 
980425 and 060218. 

4.1. Luminosity Function and Local Rate of HL-GRBs 

Following the methodology discussed above, we derive 
the luminosity function and local rate of HL-GRBs with 
our z— known HL-GRB sample. Figure 1(a) and (b) show 
the distribution of pj c _ s in (af L , a^ L )|L HL =2.3xio 52 an d 



„HL 



1 1 , ib !L )| Q HL =2 3 planes, respectively. Our sample in- 
cludes 44 bursts, and a value of p T K _ s > 0.05 marginally 
suggests the consistency between the model predictions 
and the observations. Therefore, we show the contours 
for p\_g > 0.05 only. From Fig.l we find that the pa- 
rameters are convolved. Different combinations of these 
parameters could make the consistency between both the 
luminosity and redshift distributions and the observa- 
tions. Nonetheless, strong constraints on the parameters 
could be derived from Fig.l. First, the value of af h is 
required to be smaller than 1.1 (marked by a vertical 
line in the two panels). Second, Lf h is correlated with 



a 



HL 



. This correlation is roughly described by a func- 
tion of logZ^ L /(1.2 ± 0.6)10 52 = af L3 [marked by solid 
and dashed curves in Fig. 1(a)]. Third, a£ L should be 
larger than 1.5, and a preferred value from our sample is 
2.0 < af h < 2.6 and 0.5 < af L < 0.8 [with p T K _ s > 0.15, 
marked by an ellipse in Fig. 1 (b)] . 

With these constraints on the parameters we derive 
the distribution of /5q L through simulations. We take 
«hl = o 65 ± o 15; a HL = 2.3 ± 0.3, and Z^ L /10 52 = 

(1.2 ± 0.6) x 10 Q i. We assume a normal distribution 
of the errors for 
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.,HL 



and logi^ , and simulate 
1000 sets of (af L , af L , logL^ L ). We then calculate 
pg L for each set of the parameters with Eq. 5 by tak- 
ing N = N™. The distribution of the together with 
that of these parameters are shown in Fig. 2. We ob- 
tain $ L = l.Utoil Gpc" 3 yr- 1 at a 90% confidence 
level. The 1-dimensional and 2-dimensional distributions 
of the luminosity and redshift predicted by our model with 
ai — 0.65, a 2 — 2.3, Lb — 2.25 x 10 52 erg s _1 are shown 
in Fig. 3. It is found that both the one-D z and L dis- 
tributions predicted by our model are consistent with the 
observations, and all the HL-GRBs in the 2-dimcnsional 
plane are enclosed in the 3tr region. Both $ HL and are 
consistent with those derived from the BATSE data (e.g. 
Schmidt 2001; Stern ct al. 2002; Guetta et al. 2005) and 
simulations (e.g. Dai & Zhang 2005; Daigne et al. 2006). 

4.2. LL-GRBs as the Low Luminosity End of the 
HL-GRB Population? 

We estimate p^ L with the two detections of GRB 060218 
and GRB 980425 from, 



QBepprpBepp _j_ QSwrpSv 
N = Po K<0.033 



47T 



(6) 



where Vz<o.o33 ~ 1.2 x 10 2 Gpc 3 is the volume at 
z < 0.033, n Be PP = 0.123 and fl Sw = 1.33 are the solid 
angles of the GRBM on board BeppoSAX and the BAT 
on board Swift, respectively, and T Bepp ~ 6 years and 
T Sw ~ 2 years (at the time when this paper is written) 
are the operation periods of the BeppoSAX and Swift, re- 
spectively. We obtain p^ L ~ 600 Gpc" 3 yr -1 . This is 
much larger than derived above(see also Cobb et al. 
2006 and Soderberg et al. 2006b). A simple extrapolation 
of $ HL derived above to low luminosities cannot account 
for such a high p$ h (e.g. p}f^ ~ 10 Gpc" 3 yr" 1 ; Guetta 
et al. 2004). This rate yields only a Poisson probability 
of 5.1 x 10~ 4 for detecting the two events at z < 0.0331. 
In order to argue a same population for both HL and LL 
GRBs and keep a high p}j h , the only way is to invoke a 
much larger a\ than af L derived above. We fix a 2 = 2.3, 
L b = 5 x 10 52 erg s" 1 , and N = 200, and search for a 
a\ value that can yield a local rate as ~ 600 Gpc" 3 yr" 1 . 
We obtain a.\ ~ 1.6. The derived one-D and two-D distri- 
butions of the redshift and luminosity are shown in Fig. 4. 
One can see that the model predictions are significantly 
deviated from the observations. The model greatly over- 
predicts GRBs with z ~ 1.2 and L ~ 3 x 10 52 erg s" 1 . 
Although this result could not fully rule out the possibil- 
ity that LL-GRBs are simply the low luminosity end of the 
HL-GRB population(Guetta et al. 2004), it does disfavor 
such a possibility. In the above analysis we do not consider 
the difference of the beaming effects for both HL- and LL- 
GRBs and the LF cosmological evolution effect. As we 
discuss below, the LL-GRBs are likely less collimated and 
are detectable in the nearby universe only. 

5. LL-GRBS AS A DISTINCT GRB POPULATION FROM 
HL-GRBS 

As discussed above, the high detection rate of the LL- 
GRBs motivates us to consider the LL-GRBs as a distinct 
GRB population from the HL-GRBs. The conventional 
HL-GRBs generally have a luminosity of L > 10 49 erg 
s" 1 . We therefore take a preliminary criterion of L < 10 49 
erg s _1 to select our LL-GRB sample. LL-GRBs are faint. 
They are only detectable in a small volume of the local uni- 
verse and a large amount of the population are below the 
sensitivity threshold of the detector. The observable LL- 
GRBs with Swift are rare events comparable to HL-GRBs. 
It is unlikely to establish a large sample with the current 
GRB missions, so it is difficult to investigate $ LL through 
fitting its logA^ — logP distribution or through our 1-D 
Criteria (as is done for the HL-population) . We can only 
roughly constrain the $ LL and pg L with a few detections 
and limits of LL-GRBs. GRBs 980425 and 060218 are two 
firm detections of LL-GRBs 2 . There are also two other 
marginal detections for the LL-GRBs, i.e., GRBs 031203 
(z = 0.105, L = 3.5 x 10 48 erg s" 1 ) and 020903 (z = 0.25, 
Soderberg et al. 2002; L = 8.3 x 10 48 erg s" 1 ). 

5.1. Luminosity Function and local Rate 

With the four detections and the other constraints from 
observations we constrain the LF of these LL-GRBs. The 
luminosity of these LL-GRBs range from 5 x 10 46 erg s" 1 
to 8.3 x 10 49 erg s" 1 . Assuming also a broken power law 



2 Please note that GRB 060218 shows significant hard-to-soft spectral evolution (Campana et al. 2006; Ghisellini et al. 2006) and the peak 
energy of its integrated spectrum matches the Amati- relation (Amati ct al. 2006). GRB 980425 significantly deviates this relation. Ghisellini 
ct al. (2006) argued that by considering the spectral evolution effect, GRB 980425 may be also consistent with the Amati-relation. 
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LF for the LL-population (similar to eq.[4]), we take Lf, 
around 10 47 erg s _1 and constrain ot\ and a 2 by requiring 
that the 3a contour of the 2-D distribution encloses these 
LL-GRBs. This places constraints on both ct\ and a 2 . In 
order to make the 3a contour marginally enclose the near- 
est burst, GRB 980425, but not over-predict the detection 
probability at z < 0.01, ot\ should be shallow. Similarly, 
the a 2 is constrained by GRBs 031203 and 020903. Based 
on these observational constraints we search for a\ h and 
c^ L by taking Z,£ L = (1.0 ± 0.3) x 10 47 erg s" 1 . We find 
a\ L = ± 0.5 and o^ L ~ 3.0 - 4.0 can roughly reflect 
these constraints. We use the same simulation method as 
that for HL-GRBs to derive the distribution of /3q L . The 
parameters are taken as a\ h = ± 0.5, a 2 h = 3.5 ± 0.5, 
and = (1.0 ± 0.3) x 10 47 erg s" 1 . The distribution 
of Pq L together with that of these parameters are also 
shown in Fig. 2. We obtain = 325±?^ at a 90% 
confidence level. The 2-D distribution in the (log L, log z) 
plane is shown Fig. 3. It is found that the LL-GRBs 
form a distinct "island" from the main "continental" pop- 
ulation. The detection rate of the LL-GRBs thus can be 
explained without over-predicting the HL-GRBs. These 
results suggest that the current data are consistent with 
the conjecture that LL-GRBs form a distinct population 
from HL-GRBs, with a low luminosity and a high local 
rate. The constrained luminosity functions for both HL 
and LL populations are displayed in Fig. 5 (a). 

5.2. Beaming factor 

Understanding GRB rates has profound implications for 
understanding the relation between GRBs and Type Ib/c 
SNc (e.g. Lamb et al. 2005). Compared with the local 
Ib/c SN rate (4.8 x 10 4 Gpc" 3 yr" 1 ; Marzke et al. 1998; 
Cappellaro, Evans & Turatto 1999; Folkes et al. 1999), 
the rate of LL-GRBs (on-beam only, not including those 
beamed away from the Earth) is about ~ 0.7% of the Type 
Ib/c SN rate. Most recently, Soderberg et al. (2006a) 
argued that at most ~ 10% of Type Ib/c SNe are asso- 
ciated with off-beam LL-GRBs based on their late-time 
radio observations of 68 local Type Ib/c SNc. This result, 
combined with our result, suggest that the beaming factor 
of these LL-GRBs is at most a factor of 14, in contrast 
to a higher factor (~ 100, Frail et al. 2001; Zhang et al. 
2004; Guetta et al. 2005) for HL-GRBs. This suggests 
that LL-GRBs are less collimated, with an opening angle 
typically larger than ~ 31°. This is consistent with the 
observational data of GRB 060218 (Campana et al. 2006; 
Soderberg et al. 2006b). 

6. OBSERVED LOCAL GRB EVENT RATES OF HL-GRBS 
AND LL-GRBS 

The so-called observed local GRB event rate p obs cru- 
cially depends on the luminosity function, the instrument 
threshold, and a certain redshift to estimate the rate. Fig- 
ure 5(a) shows the combined $ of both LL- and HL- GRBs 
derived from a set of ordinary parameters (a^ L = 0.65, 

af L = 2.3, L b = 1.20 x 10 52 x 10 Q ? L3 erg s" 1 ; a\ L = 0.0, 
a| L = 3.5, L\ h = 10 47 erg s" 1 ) and from two sets 
of parameters that are roughly taken as the lower limit 

(af L = 0.5, af L = 2.0, Lf L = 0.6 x 10 52 x 10"? L3 erg s" 1 ; 

a LL = _ .5, a LL = 4 0] L LL = 7x 1Q 46 crg g -l) and upper 



limit (af L = 0.8, a!? L = 2.6, Lf L = 1.8 x 10 52 x 10 Q i 
erg s" 1 ; a\ L = 0.5, c^ L = 3.0, L^ L = 3 x 10 47 erg s" 1 ) of 
the LF. It is evident that $ has two distinct components. 
We calculate the observed GRB event rates for both LL- 
and HL-GRBs as a function of "enclosing redshift" z enc 
(i.e. the volume enclosed by this redshift) for the three 
parameter sets. The results are shown in Fig. 5(b). The 
observed rates of LL-GRBs and HL-GRBs as a function 
of redshift are significantly different. p^ s keeps almost a 
constant within z < 0.1, and drops sharply at z > 0.1. At 
z ~ 0.2 p^bs i s ~ 10 Gpc~ 3 yr _1 , which is roughly con- 
sistent with that suggested by Guetta et al. (2004). The 
Pobs i s a l so almost constant within z < 0.1, but it increases 
when z > 0.1 and peaks at z ~ 2.5, being consistent with 
the observations with Swift. Although p)£ s {z) has a large 
uncertainty, (e.g. varies by almost two orders of magni- 
tude), Fig. 5(b) indicates that the detectable GRB sample 
within z < 0.4 is dominated by the LL-GRBs. 

7. DISCUSSION 

We have analyzed the luminosity function and local rate 
of HL-GRBs and extended this analysis to LL-GRBs. We 
present detailed discussion on our method, data, and anal- 
ysis results before drawing the conclusions of our analysis. 

Ideally, deriving the luminosity function requires a large, 
unbiased z— known sample. The current GRB sample, 
however, is still small and suffers from some observational 
biases. First, the trigger probability of a burst depends on 
its peak flux. A weaker burst near the instrument thresh- 
old has a lower trigger probability. This effect affects the 
completeness of the sample at a given instrumental thresh- 
old. As shown in Fig. 3, a gap between the data points 
and the threshold curve is observed. The gap manifests 
this bias. Second, the redshift measurements also have bi- 
ases against weak GRBs (e.g. Bloom 2003). In order to 
evaluate whether these observational biases could signifi- 
cantly affect our results, we use the 2-D distribution in the 
(log L, log z) plane at 3a significance level to examine our 
results. If these biases do not ruinously affect our results, 
the 3cr contour of the GRB distribution in the (log L, log z) 
plane should enclose the current GRB sample; otherwise, 
the model should predict many undetected GRBs near the 
threshold and the 3cr contour cannot cover the scattered 
data point in the (log L, log z) plane. From Fig. 3 we find 
the data points are enclosed by the 3a contour for the HL- 
GRBs, suggesting that these observational biases do not 
significantly affect our results for HL-GRBs. However, the 
current Swift GRB sample with redshift measurements is 
still small. One prominent issue is whether both the red- 
shift and luminosity distributions derived from the cur- 
rent Swift GRB sample represent the true observed ones 
for the threshold of Swift. In order to address this issue 
we plot also the GRBs detected by previous GRB mis- 
sions in Fig. 3. It is found that most of these GRBs are 
enclosed by the 3a contour defined with the current Swift 
HL-GRB sample. Two GRBs, 030329 and 040701, are out 
of the 3a contour. In order to make a 3a contour to en- 
close the two GRBs, we find that the luminosity function 
should have ot\ > 1.05, a 2 ~ 3.0, and Lj, ~ 6 x 10 52 erg 
s _1 . We also show that the 3a contour and the predicted 
L and z distributions of this luminosity function in Fig. 3 
(dashed line). It predicts more GRBs with L ~ 10 49 - 10 51 
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erg s^ 1 (roughly the low luminosity end of the HL-GRBs) 
than observations by Swift. If this is indeed an observa- 
tional bias for the HL-GRB sample with Swift, the true oti 
should be larger than that derived from the current Swift 
GRB sample, say, «i - 1. With such a LF, GRB 020903 
is enclosed in the 3cr contour and 031203 also marginally 
belongs to the HL-GRB group. However, even with such 
a LF, GRBs 980425 and 060218 are definitely not merged 
into the HL-GRB group. 

Another effective approach to evaluate the biases is to 
fit the observed log A — logP distribution. The current 
Swift sample has only ~ 200 bursts. It still cannot be 
used for making such a fitting. A sample of <~ 500 bursts 
with ~ 120 rcdshift measurements could be obtained by 
Swift in a period of 5 years. Such a sample could refine 
the constraints on the LF of HL-GRBs by incorporating 
our approach with the log N — logP fitting 3 . 

A cosmological evolution effect may influence the deriva- 
tion of the luminosity function (e.g. Lloy-Ronning et al. 
2002). So far there is no robust evidence of such an evo- 
lutionary effect with the GRB sample with spectroscopic 
redshifts. This effect also does not affect $ LL , since the 
LL-GRBs are likely detectable only at low redshifts. It 
may influence the derivation of <1> HL . A large sample of 
spectroscopic redshift measurements is required to more 
robustly evaluate this effect. The current HL-GRB sample 
is inadequate for this task. Therefore, we do not consider 
this effect in this analysis. In the future, using a sam- 
ple established by Swift within 5 year operation one may 
be able to robustly model this effect, and the constraints 
on the luminosity function of HL-GRBs could be further 
refined. 

Our results tentatively suggest that LL-GRBs are an 
intrinsically different GRB population, characterized by a 
high local rate, low luminosity, and small beaming factor 
comparing to the HL-GRBs. The detectable low redshift 
GRBs should be dominated by the LL-GRBs. Le & Der- 
mer (2007) draw the similar conclusion through modeling 
the redshift distribution of the Swift GRBs. However, we 
pointed out that the current data do not conclusively ex- 
clude the possibility that only one population of GRBs 
is responsible for both HL- and LL-GRBs(Guetta et al. 
2004). 

There are two scenarios to explain the nature of the 
LL-GRBs. One scenario is that these GRBs are standard 
HL-GRBs viewed off-axially (e.g. Nakamura et al. 1999; 
Yamazaki et al. 2003). In order to account for the step- 
like two-component LF displayed in Fig. 4, the jet must 
include two distinct components, i.e. a narrow HL compo- 
nent and a very wide LL component. Such a jet configura- 
tion is different from the conventional jet-cocoon picture in 
the standard collapsar model in which the cocoon compo- 
nent is not as broad as 31° (e.g. Zhang et al. 2003). The 
lack of detection of radio re-brightening in GRB 980425 
and GRB 031203 (Soderberg et al. 2004a,b) has also 
greatly constrained this scenario (Waxman 2004a, b). 

The second scenario is that LL-GRBs are intrinsically 
different from HL-GRBs (e.g. Wang et al. 2000; Soder- 
berg et al. 2004b; Mazzali et al. 2006). Actually, the 
LL-GRBs as a unique GRB population has been proposed 

3 The log N — logP distribution is sensitive to the instrument thre 
such a fit. CG-RO/BATSE established a large sample and such a fit gi 
with our results shown in Fig.l 



since the discovery of GRB 980425/SN1998bw with its ob- 
served under-luminous, long-lag, simple light curve fea- 
tures (e.g., Soderberg et al. 2004b). Stern et al. (1999) 
suggested a group of "simple" bursts with peak fluxes near 
the BATSE trigger threshold. Norris (2002) found that 
the proportion of long-lag, single-pulse bursts (similar to 
the three LL-GRBs) within long-duration bursts increases 
from negligible among bright BATSE bursts to^ 50% at 
the trigger threshold, and their peak flux ~ 2 orders of 
magnitude lower than that of the brightest bursts. Taken 
together with the fact that three nearby GRBs, 980425, 
031203, and 060218, are long-lag and under-luminous (e.g., 
Sazonov et al. 2004; Liang et al. 2007), an intuitive specu- 
lation of the observed LL-GRBs is that that they are prob- 
ably relatively nearby (e.g., Norris et al. 2005) and their 
local rate should be much higher that conventional HL- 
GRBs as we derived in this analysis. This scenario calls for 
a different type of progenitors (e.g. neutron stars rather 
than black holes) for LL-GRBs from those of HL-GRBs 
(e.g. Mazzali et al. 2006; Soderberg et al. 2006b). The 
abundant multi- wavelength observations of GRB 060218 
lead us make a step toward to the nature of this kind of 
events (Wang & Meszaros 2006; Fan et al. 2006; Dai et 
al. 2006; Li 2007; Ghisellini et al. 2006,2007; Wang et al. 
2006; Gupta & Zhang 2007; Murase et al. 2006; Amati et 
al. 2006; Liang et al. 2007). Future detectors more sen- 
sitive to Swift BAT, e.g. EXIST (Grindley 2006), would 
detect much more of these events and greatly increase the 
sample of LL-GRBs. 

8. CONCLUSIONS 

We have derived the luminosity function and the lo- 
cal rate of the HL-GRBs with the Swift z-known GRB 
sample and explored the possibility that LL-GRBs as a 
unique GRB population from the HL-GRBs. We find that 
Pg L is ~ 325i 3 y7 Gpc~ 3 yr _1 , which is much higher than 

p$ h (1.12±%il Gpc~ 3 yr- 1 ). This rate is - 0.7% of the lo- 
cal Type Ib/c SNe. Our results, together with the finding 
that less than 10% of Type Ib/c SNe are associated with 
off-beam GRBs, suggest that LL-GRBs have a beaming 
factor typically less than 14, or a jet angle typically wider 
than 31°. The high local GRB rate, the small beaming fac- 
tor, and low luminosity make the LL-GRBs distinct from 
the HL-GRBs. Although the current data could not fully 
rule out the possibility that both HL- and LL-GRBs are 
the same population, our results suggest that LL-GRBs 
are likely a unique GRB population and the observed low 
rcdshift GRB sample is dominated by LL-GRBs. 

We note that Guetta & Delia Valle (2007) reproduce the 
three main results, i.e., the high local rate of LL-GRBs, 
HL- and LL-GRBs as two populations, and wide jet open- 
ing angle for LL-GRBs, three months after we posted this 
paper to astro-ph. They also suggest that the possibility 
of one population for both the HL- and LL-GRBs (Guetta 
et al. 2004) cannot be convincingly ruled out. 

We thank Don Lamb for helpful discussion. This work 
was supported by NASA under grants NNG05GB67G, 
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>ld. One has to use a sample observed by a given instrument to make 
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Fig. l. — Distributions of Pk-s m the planes of (a^ L , ^b IL )laH L =2.3 (panel a) and (ctf h , )|j,hl_ 2 3x10 52 (panel b) 
derived with the current Swift HL-GRB sample. The vertical dashed line marks a limit on oi\. The curves in the panel 
(a) show the relation and its errors between Lj^ L and . The ellipse in panel (b) marks the best parameters of af L 
and af L . 
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Fig. 2. — Distributions of the GRB local event rates and the parameters of the LFs for LL-GRBs (left panels) and 
HL-GRBs (right panels). 
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Fig. 3. — Jointed contours of the logarithmic GRB detection rate [log(dN / dt)] distribution in a 2-dimensional (2-d) 
[logL, log(z)]-plane as compared with observational data [panel (a)], assuming that the HL- and LL- GRBs are two 
distinct populations. The two firm LL-GRBs are denoted by stars, and the Swift HL-GRBs are denoted by solid dots. 
The hatched region marks the limitation of the Swift/BAT detectability, where the threshold is derived by using the 
Swift/BAT sensitivity in the 50-150 keV band for a standard GRB with E p = 200 keV in the GRB local frame. The 
bold-faced, solid curves in panel (a) marks the 3<r confidence level of the 2-d distributions for the HL- and LL-GRBs. The 
comparisons of the observed 1-d distributions of log L and log z with the model predictions are presented in panels (b) and 
(c), respectively. The dashed curve in the panel (a) and the dashed lines in the panels (b) and (c) are, respectively, the 
3cr contour of the 2-d distribution and the corresponding 1-d distributions derived from a LF with = 1.05, aJ? L = 3, 
and L^ L = 6 x 10 52 erg s _1 , which gives a 3cr contour that can enclose all the HL-GRBs observed by Swift and pie-Swift 
missions (see §7 in the text). 
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Fig. 4. — Same as Fig. 3 but for the case that the HL- and LL- GRBs are assumed to belong to the same population. 



11 




Fig. 5. — Panel (a): The combined LFs of both LL- and HL- GRBs derived from a set of ordinary parameters (solid line) 
and from two sets of parameters that are roughly regarded as the lower (dash-dotted line) and upper (dashed line) limits 
of the LFs. Panel (b): The observed GRB event rates for both LL- and HL-GRBs as a function of "enclosing redshift" 
z enc (i.e. the volume enclosed by this redshift) for the three parameter sets shown in panel (a). Same line styles for 
different models are adopted in both panels. 



